Nanoscale Phase Coexistence and Percolative Quantum Transport 



Sanjeev Kumar and Pinaki Majumdar 
Harish- Chandra Research Institute, 
Chhatnag Road, Jhusi, Allahabad 211 019, India 
(Nov 15, 2003) 



^1" 

o 
o 

(N 
m 



We study the nanoscale phase coexistence of ferromagnetic metallic (FMM) and antiferromagnetic 
insulating (AFI) regions by including the effect of AF superexchange and weak disorder in the double 
exchange model. We use a new Monte Carlo technique, mapping on the disordered spin-fermion 
problem to an effective short range spin model, with self-consistently computed exchange constants. 
We recover 'cluster coexistence' as seen earlier in exact simulation of small systems. The much larger 
sizes, ~ 32 x 32, accessible with our technique, allows us to study the cluster pattern for varying 
electron density, disorder, and temperature. We track the magnetic structure, obtain the density of 
states, with its 'pseudogap' features, and, for the first time, provide a fully microscopic estimate of 
the resistivity in a phase coexistence regime, comparing it with the 'percolation' scenario. 
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The issue of multiphase coexistence in transition metal 
oxides has been brought to the fore by a set of remark- 
able recent experiments on the manganites [1,2]. These 
experiments probe the atomic scale magnetic correla- 
tions [3], structural features [4], or 'conducting' prop- 
erties, i.e tunneling density of states [5], through local 
spectroscopy. The bulk thermodynamic and transport 
properties of these systems, including the 'colossal mag- 
netoresistance' (CMR) seem to have a correlation with 
the physics of electronic phase coexistence as visualised 
in these nanoscale experiments. 

Establishing a first principles theoretical connection 
between cluster coexistence and 'anomalous' transport, 
however, has been difficult. Only real space Monte Carlo 
(MC) techniques [6] allow access to the cluster phase, and 
accessible sizes, ~ 8 x 8 in two dimension (2d), are too 
small to study transport. The approach in current use 
[7] is phcnomcnological, using classical resistor networks 
to model the data within a 'percolation scenario'. While 
'percolation' surely plays an important role in transport, 
standard percolation theory [8] is inadequate [9] to ex- 
plain the data. More fundamentally, such an approach 
does not explain the temperature and field dependent re- 
sistance of the 'building blocks', or correlations and hys- 
teresis in the cluster distribution, or the dependence of 
bulk resistance on 'spin overlap' between clusters. 

We address these questions within a 'spin-fermion' 
model in this paper, employing a new MC technique [10] 
that allows access to system sizes ~ 32 x 32 in 2d. This 
can probe the regime L$ 3> L c 3> ao, where Ls, L c 
and ao are system size, typical cluster size, and lattice 
spacing, respectively. We study the coexistence regime, 
obtain the typical cluster size, and calculate the spec- 
tral density and conductivity in the mixed phase. To our 
knowledge this is the first microscopic calculation to clar- 
ify the connection between nanoscale phase coexistence 
and transport in a fully quantum mechanical itinerant 
electron system. 



The coexistence of two phases with distinct electronic, 
magnetic, and possibly structural, properties is best con- 
ceived at T — in a clean system. For a specified 
chemical potential, fj,, the ground state configuration of 
the spin and lattice variables, {Xi} say, assumed clas- 
sical, is determined by SS/SXi = 0, where £{A^;/i} is 
the energy of the system in the {Xi} background. The 
minimum, £ m in(h i ), usually occurs for a unique {Xi} at 
each jj,, and in this background the electron density n(/j) 
is also unique. However, in the presence of competing 
interactions, two distinct {Xi} configurations could be 
degenerate minima of E at some /i. The corresponding 
\i = fx c marks a first order phase boundary, and the two 
'endpoint' densities, n\ and n 2 bracket a region of co- 
existence. There is no homogeneous phase with density 
between n\ and n 2 . In this regime the system breaks up 
into two macroscopic domains, with density rv\ and n 2 . 

It was pointed out by Dagotto and coworkers [6], and 
known earlier in the context of classical spin models [11], 
that disorder is a singular perturbation in these sys- 
tems. Even weak disorder breaks up the 'macroscopic' 
domains into interspersed locally ordered clusters of the 
two phases. This nano/mesoscopically inhomogeneous 
strongly correlated phase can be studied only within a 
real space framework. 

In this paper we study the key qualitative issues of 
phase coexistence in itinerant fermion models by consid- 
ering the competing effects of double exchange (DE) and 
superexchange (SE) in the presence of weak disorder. We 
consider H — H e i + Haf in 2d, with: 



H, 



u 



(e» - (j,)rii 



Jh S,.(Tj 



(1) 



and H A f = JsJ2(ij) S i- S j- Tne hopping t i3 = -t for 
nearest neighbours, and ti is uniformly distributed be- 
tween ±A/2. We assume Jjj/t — > oo. The parameters 
in the problem are A/t, Js/t and density n (or chemical 
potential /i). We use |S»| = 1. For Jh /t — > oo, the Hamil- 
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tonian, in the projected basis [10], assumes a simpler 
form: H ei = -*E<ij) A?( ^ihi + h - c ) + E*( e * ~mH- 
The hopping amplitude, = fije l ® ij , between locally 
aligned states, can be written in terms of the polar angle 
(6i) and azimuthal angle of the spin Si as, cos ^-cos^- 
+sin^sin^fe~ l ^ i_< ^). The 'magnitude' of the overlap 
is fij = + Sj.Sj)/2, while the phase is specified by 
tan$ij = Im(gij)/Re(gij). 

To make progress we need the 'effective Hamiltonian' 
controlling the Boltzmann weight for the spins. For- 
mally, this is: H eff {Si} = -(1/ 'f3)logTre- p{H ' l+HAF) . 
The 'exact' MC directly generates equilibrium spin con- 
figurations of H e f f through iterative diagonalisation [6] . 
We make the approximation: H e ff = — zC(ij) ^ijfij + 
Jg Y]/jj\ Sj.S,-, with D i3 being determined self consis- 
tently [10] as the average of = (e l ^ iir )\^lj + h.c) over 
the assumed equilibrium distribution. This approxima- 
tion has been extensively benchmarked for the clean DE 
model [10] and we will put up similar comparisons for the 
present problem in the near future [12]. 

Although the Dij enter as 'nearest neighbour' ex- 
change, they arise from a solution of the full quantum 
statistical problem in the disordered finite temperature 
system. In the presence of competing interactions and 
quenched disorder, this leads to a set of strongly inho- 
mogeneous, spatially correlated, and temperature depen- 
dent 'exchange' D^. 

At consistency, fermionic properties are computed and 
averaged over equilibrium spin configurations. We work 
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FIG. 1. Magnetisation, m(T), and structure factor, 5* q (T) 
at q = {-k,tt}. Panel (a) — (b): data for A = 0.2, panel 
(c) — (d): data for A = 1.0. The legends are common to 
panels (a) — (d). Panel (e): The characteristic temperature 
(see text), T ch (n, A): triangle down: AF at A = 1.0, square: 
AF at A = 0.2, triangle up: FM at A = 1.0, circles: FM 
at A = 0.2. Panel (/): m(T = 0) with varying n, circles: 
A = 0.2, squares: A = 1.0. System size: 24 x 24. 




FIG. 2. Cluster pattern: Rows 1-2 arc for A = 0.2, row 
3-4 are for A = 1.0. Mean density n = 0.1 in both cases. 
Upper row in each set shows the thermally averaged density 
profile, n r , the lower row shows the nearest neighbour spin 
correlation: /| = (S(r).S(r+<5)). Left to right along each row, 
the T variation is T/T c ~ 2.0, 1.0, 0.0. For n r dark regions are 
high density. For /f dark regions are FM, white regions AF. 
System size 32 x 32, data for a specific {e^} in each case. 

at specified mean density, fixing /j, through iteration. 
Transport properties are computed using the Kubo for- 
mula [10], employing sizes <~ 24 x 24 to 32 x 32, and 
averaged typically over 30 realisations of disorder, with 
averaging over ~ 50 equilibrium spin configurations at 
each T for each realisation of {e^}. Since the d.c con- 
ductivity is not directly accessible in a finite system, we 
compute the finite frequency conductivity at the scale 
of mean level spacing, uol cx 8t/L 2 , with ul — 0.04t at 
L = 32. The conductivity results are in units of (ne 2 )/Ji. 
We have checked the adequacy of disorder average and 
stability with respect to system size variation [12]. 

Our principal results include the direct visual evidence 
of cluster coexistence (Fig. 2), the disorder dependence of 
spectral density (Fig. 3) and the resistivity (Fig. 4) in the 
coexistence regime. 

At low electron density, the competetion in the DE+SE 
model is between a {n, n} AF phase and a ferromagnct 
(FM). We set Js = 0.05 and scanned in \x to locate the \i c 
for the first order boundary. The density changes from 
n = to n <~ 0.20 at the discontinuity. These results 
were obtained within the H e f / scheme and cross checked 
with exact MC. Moderate disorder smooths out the dis- 
continuity in n(/i) converting it to a sharp crossover [12]. 
In this paper we fix Js — 0.05 and probe the coexistence 
regime, n ~ — 0.2, with varying electron density, tem- 
perature (T), and A = 0.2 and A = 1.0. All the results 
are obtained by cooling the system from the high T 
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FIG. 3. Low energy density of states, at n = 0.10. Panel 
(a): data at A = 0.2, panel (6): data at A = 1.0. The 
DOS is plotted with respect to w - p(T), incorporating the 
T dependent shift in p. Lorentzian broadening of 8 functions 
to 0.03. Data averaged over 30 copies. 

paramagnetic, approximately homogeneous, phase. 

Fig.l shows the magnetisation, m(T), and the AF peak 
in the magnetic structure factor ,S q at q = {ir, w} to il- 
lustrate the evolution from the AF to the FM phase with 
increasing n. Panel (a) — (b) are at A = 0.2 and panel 
(c) — (d) at A = 1.0. The 'extremal' densities are either 
strongly AF or FM, while for n ~ 0.08-0.12 both FM and 
AF reflections have finite weight. Panel (e) tracks the 
'characteristic temperature', T c h(n, A), identified from 
the maximum in d 2 OjdT 2 where O is the appropriate 
order parameter (of the FM or AF phase) . In the shaded 
region, 0.10 < n < 0.14, it is difficult to resolve the T c h 
accurately. Panel (/) shows the change in 'saturation 
magnetisation' m(T = 0) with increasing n and chang- 
ing disorder. The small moment regime, m < 0.1, for 
n < 0.1 is a 'ferro- insulator' phase, as we will discover 
from the transport data. In this regime the moments in 
different clusters are only weakly correlated. 

Fig. 2 provides direct visual evidence of 'clustering', the 
top two panels showing thermally averaged local den- 
sity and nearest neighbour spin correlation for T/T c ~ 
2.0,1.0,0.0 (left to right), at A = 0.2. The lower panels 
show data at A = 1.0. We will quantify the 'typical scale' 
of clusters later, but from the density distribution it is 
obvious that the pattern for A = 1.0 is more fragmented 
than for A = 0.2. The 'density contrast' reduces with 
increasing T, as the spins in the AF regions fluctuate out 
of antiparallcl alignment, and the carriers can partially 
explore these regions. The spin correlations also weaken, 
with FM (black) and AF (white) regions giving way to 
an uncorrelated (grey) background. 

Fig. 3 shows the low energy density of states (DOS) for 
n = 0.1. Even at weak disorder, panel (a), there is a 
'pseudogap' in the system averaged DOS, at the lowest 
temperature, T ~ 0.01. However, the pseudogap rapidly 
fills up with increasing T, and the DOS tends towards 
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FIG. 4. Resistivity, p(T), with varying n. Panel (a): 
A = 0.2, panel (&): A = 1.0. Results obtained on 'cooling'. 
System size 24 x 24, data obtained by inverting the 'mean 
conductivity'. Average over 20 — 50 copies of {e^} (error bars 
comparable to symbol size, unless otherwise indicated). Re- 
sults on size 32 x 32 are similar. Inset, panel (b): Normalised 
T = conductivity, at A = 1, -vs- FMM surface area. 



the universal profile of the spin disordered 2d DE model. 
At stronger disorder, panel (b), the 'dip' in the DOS at p, 
is deeper. Due to stronger pinning, the clusters, and the 
gap feature, survives to higher T, and only for T 3> T c h 
tends to the asymptotic form. While this data focuses on 
the system averaged DOS, which can be probed by pho- 
toemission (PES), tunneling spectroscopy would track 
the local DOS, which shows [12] a clean 'gap' in the AF 
regions, and finite DOS in the (larger) metallic clusters. 

Fig. 4 shows resistivity, p(T), with varying n and dis- 
order. The comparison of panel (a), A = 0.2, and panel 
(6), A = 1.0, indicates that p(T) for stronger disorder is 
systematically larger. The trends, however, are similar 
in the two cases and allows a tentative classification of 
the 'global' aspects of transport. At both A = 0.2 and 
A = 1.0 there is a critical density, n c (A), below which 
p(T = 0) diverges, indicating the absence of any con- 
nected 'conducting path'. The inset to panel (b) shows 
the trend in the T = conductivity with increasing FM 
cluster area, ofm- 

We suggest the following qualitative picture of trans- 
port in the coexistence regime based on the data in Fig. 4. 
(a). At low T, and for n > n c (A), there is, by definition, 
some connected 'metallic path' through the sample. It 
is reasonable to assume that at least in this connected 
region the spins are aligned due to DE. Since cluster di- 
mension and electron wavelength are comparable in our 
system a large fraction of the resistance arises from the 
non trivial geometry of the current path. This scattering 
is quantum mechanical, and the low T metallic phase cor- 
responds to a quantum percolative regime [13], for roughly 
0.15 > n > 0.08 at A = 1.0. (6). With increasing T, the 
resistance of the conducting network increases due to DE 
spin fluctuations and, till the network is disrupted, there 
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FIG. 5. Typical FMM cluster size: (a) A = 0.2 and varying 
T, (b) A = 1.0, varying T. The inset to panel (6) highlights 
the reduction in cluster size with increasing disorder. 

is a regime dp/dT > 0. This 'metallic' behaviour oc- 
curs despite the very large residual resistivity. This is 
a regime of weak magnetic scattering on the conducting 
network. If n 3> n c , e.g, n — 0.20 in Fig. 4. (b), so that 
inhomogeneities are weak, p(T) will smoothly increase to 
the T T c h asymptotic value, (c). With further increase 
in T, in the n>n c regime, the spin disorder can destroy 
some of the 'weak links', disrupting the conducting net- 
work and leading to a sharp increase in p(T) see, e.g, 
n = 0.1 in Fig. 4. (b). This correlates well with reduction 
in typical cluster size with increasing T, discussed in the 
next paragraph. Depending on n and the extent of dis- 
order (and system size, in a simulation) there could be 
a rapid rise or a 'first order' metal-insulator transition 
(MIT). This is spin disorder induced breakup of clusters 
driving a MIT. (d). Beyond this 'MIT' the conduction 
is through the 'insulating' regions, with isolated patches 
contributing to nominally activated transport, as visible 
in Fig. 4. (b) for n < 0.12. For T > T ch , as the structures 
disappear, Fig. 2, and the system becomes homogeneous, 
p(T) is dominated by spin disorder scattering. This is a 
diffusive regime with saturated spin disorder scattering, 
visible for T > 0.08 at all densities, (e). For n < n c , the 
regime of low density isolated clusters, p(T) falls mono- 
tonically, and the response is typical of low density ferro- 
magnetic polarons in a AF background. This occurs for 
n < 0.08 at A = 1.0, and at lower n at A = 0.20. 

Finally, Fig. 5 shows the typical size of FMM clusters, 
inferred from a Lorentzian fit to the magnetic structure 
factor, i.e, 5 q ~ (q 2 +£~ 2 ) -1 . The resulting correlation 
length depends on n, A and T, decreasing with increasing 
disorder and T, and increasing with increasing density. 
The main panels, (a) and (6) , highlight the n dependence 
at different T and A, while the inset in Fig.5.(b) replots 
the same data to highlight the dependence on disorder. 
The dependences, overall, are intuitive, and strengthen 
the proposed transport scenario. 

Our results on coexistence share several generic fea- 



tures with the manganites [3-5], but there are key dif- 
ferences too. Apart from the difference between 2d and 
3d, these are (i) the important regime of coexistence in 
manganites is between FMM and charge ordered insu- 
lator (COI), and crucially involves the Jahn- Teller (JT) 
phonons. This requires enlarging our model, (ii) It is 
supposed that the FMM and COI regions are nominally 
of equal density, which is why such domains can survive 
over pm scale. Our clusters are 'charged' due to the den- 
sity difference between FMM and AFI. In any real system 
Coulomb effects [14] would have to be considered in such 
a case. We will discuss such effects separately [12], as well 
as the effect of an applied magnetic field on transport. 

In conclusion, we have presented microscopic results 
on quantum transport across a regime of phase coexis- 
tence in an itinerant fermion model, and correlated it 
with thermal evolution of the spatial structures. The 
method readily generalises to the two orbital model with 
JT phonons, to be discussed in the near future. 
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